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Abstract 

We demonstrate the growth of graphene nanocrystals by molecular beam methods that employ 
a solid carbon source, and that can be used on a diverse class of large area dielectric substrates. 
Characterization by Raman and Near Edge X-ray Absorption Fine Structure spectroscopies reveal a 
sp 2 hybridized hexagonal carbon lattice in the nanocrystals. Lower growth rates favor the formation 
of higher quality, larger size multi-layer graphene crystallites on all investigated substrates. The 
surface morphology is determined by the roughness of the underlying substrate and graphitic 
monolayer steps are observed by ambient scanning tunneling microscopy. 

PACS numbers: 81.10.Pq,61.46.Hk,78.35.+c,61.05.C- 
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I. INTRODUCTION 



Realization of the enormous potential of graphene in condensed matter science and nan- 
otechnology requires use of state-of-the-art processing methods. While mechanical exfo- 
liation of graphite offers ultra-high quality graphene flakes on arbitrary substrates, it is 
not suitable for creation of large-area films required for many experimens and for poten- 
tial applications ~ Thermal decomposition of SiC and chemical vapor deposition (CVD) on 
transition metals result in large area material, but the methods are limited to specialized 
substrates.^ 

Molecular-beam epitaxy (MBE) enables growth of precisely tailored semiconductor struc- 
tures that have greatly advanced fundamental condensed matter science.^ MBE growth of 
superb semiconductor quantum structures is facilitated by access to suitable substrates and 
by availability of atomic or molecular beams from ultra-clean sources in an ultra-high vac- 
uum (UHV) environment. The lack of a substrate with suitable interface properties for 
epitaxy hinders successful MBE-growth of graphene. MBE-like growth of graphene layers 
on arbitrary substrates would enable electrical and magnetic doping and the creation of 
sophisticated graphene-based heterostructures for fundamental research. Consequently, de- 
velopment of MBE growth of graphene on various substrates is currently an area of intense 
research activityP^ 

We report the fabrication of large area conducting graphene films by a MBE-inspired 
method that we call molecular-beam-growth (MBG). In MBG the films are directly grown 
on the substrate surface by a carbon beam that is created from a solid carbon source in an 
UHV environment. We demonstrate that MBG films of graphene nanocrystals can be grown 
on a diverse group of dielectric substrates. The combination of highly controllable growth 
conditions and dielectric substrates produces films that do not require exfoliation for further 
experiments, and facilitates comprehensive in-depth characterization. 

We employ Raman scattering spectroscopy, Near edge X-ray absorption fine structure 
spectroscopy (NEXAFS), and scanning tunneling microscopy (STM) to explore the key 
parameters required to obtain good quality ultra-thin graphitic films by MBG. Raman spec- 
troscopy probes the crystallinity and overall quality of the layers, and shows that the carbon 
atoms are arranged in a hexagonal lattice. NEXAFS spectroscopy is used to characterize the 
nature and geometry of the carbon bonds and overall quality of the films, including sp 2 :sp 3 
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ratios.^ The NEXAFS data unequivocally demonstrate planar layered sp 2 graphitic bonds 
in films grown under optimized conditions. The STM measurements reveal monoatomic 
steps and electrical conductivity. 

The quality and size of the graphene nanocrystals in the MBG films strongly depend 
on the growth conditions. Most notably, we find that the growth rate (GR), in particular, 
plays a crucial role. If the GR is too high the sp 2 carbon bonds in the nanocrystals form a 
three-dimensional (3D) arrangement that precludes the layered structure of graphene multi- 
layers. By lowering the GR (below 0.5A/min) we successfully grow layered MBG films that 
show clear signatures of a two-dimensional lattice of sp 2 -carbon in NEXAFS and Raman 
spectroscopy. On the basis of these results we discuss the MBG growth mode of ultra-thin 
films of graphene nanocrystals. 

The paper is organized in sections. Section II. A describes growth methods and the con- 
ditions for growth on diverse substrates. Section II. B describes the experimental analytical 
tools used to probe the ultra-thin MBG films. The experimental results and their evalua- 
tion are presented in Part III: Section III. A presents Raman and NEXAFS results of MBG 
growth on Si02 substrates. Section III.B describes NEXAFS experiments obtained on var- 
ious substrates. Section III.C presents NEXAFS results from thick MBG films. The last 
section, III.D, describes ambient STM measurements of the conducting surface of the MBG 
ultra-thin graphene films. Estimates, based on Raman scattering data, of the sizes of the 
graphene nanocrystals are presented in Part IV. Part V concludes this paper. 

II. EXPERIMENTAL METHODS 
A. Molecular beam growth 

MBG film growth of graphene nanocrystals on dielectric substrates, such as amorphous 
Si02 (300 nm), crystalline Mica and crystalline Si, is achieved in the set-up illustrated in Fig. 
[T] The substrates are cleaned by sonication in acetone and isopropanol prior to loading in the 
growth chamber. Mica samples are cleaved ex-situ and are loaded immediately into the UHV 
system. The UHV chamber (base pressure ~ 6 x 10~ 10 mbar) incorporates a solid carbon 
source that is made of glassy carbon, similar to that employed for p-doping in III-V MBE 
and other carbon-related growthJ^HIE] The dimensions of the carbon source are 10x2.5x0.3 
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mm 3 . The carbon source is heated by a DC current of ~15 A to an operating temperature 
of ~2100°C which is monitored by a Marathon MM Raytech optical pyrometer. The solid 
carbon source is located in close proximity to the substrate, as shown in Figs. [IJa) and (c). 
The substrates are heated to ~ 400°C to remove adsorbed water before the growth. The 
maximum pressure reached during growth is ~ 5 x 1CT 8 mbar. Due to the proximity of the 
solid carbon source, the temperature of the substrates during growth reaches ~ 500°C, with 
a gradient of less than 100°C over their 25 mm length. 

In the growth set-up shown in Fig. 1(a), D is the distance between the carbon source 
and the sample (~15 mm) and d is the position on the substrate. Go is the thickness at 
d = 0. In this configuration the flux of carbon atoms is relatively high at the near end of 
the substrate (d = 0) and decreases significantly along the length of the substrate. The 
geometrical dependence of the flux is best described as a GR gradient along the length 
of the substrate.^ The calibration of the GR is achieved by measuring the profile of a 
thick MBG films (> 30 nm) on a Si02 substrate using an atomic force microscope or optical 
profilometer. The position-dependent GR(d), derived from the position-dependent thickness 
0(d), is calculated according to the expression^ 

where t is the deposition time. The maximum GR (typically 1-2 A/min) is reached for 
d = 0. As d increases GR(g() decreases to a minimum value of 0.1 A/min, or less. 



B. Characterization methods 



Raman spectroscopy is a widely used technique for characterization of carbon-based 
materials.^231 NEXAFS provides a direct, element-specific probe of bond type and ori- 
entation with a high surface sensitivity that enables evaluation of sp 2 :sp 3 -bond ratios and 
the degree of planarity of ultra-thin (single layer) films. 15 Since sp 2 -hybridized carbon layers 
have unique spectral fingerprints in both Raman and NEXAFS spectroscopies, the combi- 
nation of these two methods is particularly suited to probing the crystallinity, bond type 
and bond configurations (2D vs 3D) of the ultra-thin MBG films. 

For the Raman experiments a Renishaw in Via micro- Raman set-up, equipped with a 
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movable x-y-z stage was employed. The laser power was set to less than 3 mW and was 
focused with a 100 x lens to a spotsize of ~ 0.5/im. 

Carbon Is NEXAFS measurements have been performed at the NIST beamline U7A of 
the National Synchrotron Light Source (NSLS). Measurements were performed in partial 
electron yield (PEY) mode with a grid bias of -200 V, selected to optimize the surface sensi- 
tivity of the measurement and thereby the signal from the graphene film. Angle-dependent 
NEXAFS was obtained by changing the angle between the incoming x-ray beam (and there- 
fore the E-field vector) and the sample between 20°and 70°, corresponding roughly to out- 
of-plane and in-plane bond resonances, respectively. The reference absorption intensity (Io) 
of the incoming x-ray beam, measured on a gold coated mesh positioned just after the refo- 
cusing optics, was measured simultaneously and used to normalize the spectra to avoid any 
artifacts due to beam instability. A linear background was subtracted from a region before 
the absorption edge (278-282 eV). Spectra were normalized by area with respect to carbon 
concentration using a two-point normalization: area normalization between 282 and 300 eV 
and a continuum normalization in the region 330-335 eV (atomic normalization). 

Ambient STM and atomic force microscopy (AFM), in tapping mode, have been per- 
formed to get additional insight into the surface morpholgy of the grown films. 



III. EXPERIMENTAL RESULTS AND DISCUSSION 



Figure [2] shows typical Raman and NEXAFS measurements for a MBG film grown on 
a 300 nm-thick Si0 2 layer on Si. As discussed below, these results show, unambiguously, 
the formation of graphene nanocrystals along the whole gradient and reveal that the grown 
layers have different structural properties depending on the GR. Two main regions with 
very distinct characteristics can be identified. The dashed line in Fig. [2] marks the border 
between these two regions, corresponding to high GR (upper half) and low GR (lower half), 
respectively. 
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A. MBG films on amorphous SiC>2 substrates 

1. Raman spectroscopy 

Characteristic Raman signatures of optical phonons for graphite are observed along the 
GR gradient, as displayed in the color plot of Fig. [2] (b).^^3 The band at ~ 1600 cm -1 
results from superposition of the G and the D' modes. The G mode is a long wavelength 
optical phonon originating from in-plane bond-stretching motion of pairs of sp 2 hybridized 
carbon atoms. The D' mode is induced by disorder and requires intra-valley electron-phonon 
scattering. The D mode at ~1344 cm -1 , which requires the presence of six-fold aromatic 
rings, is induced by disorder, such as edges or atomic defects. The bands resolved at higher 
Raman shifts are also well known: The 2D mode at ~2700 cm -1 (a.k.a. G'), the G+D band 
slightly below ~3000 cm -1 and a third one at ~3200 cm -1 matching the energy of the G+D' 
mode. We observe the latter, but it is not discernible in the color plot of Fig. [2] (b). 

The intensity of all Raman features decreases with decreasing GR (film thickness), while 
the relative intensities of D and G bands vary with the GR: For higher GR (upper part 
of Fig. [2]), the D mode is more intense than the G mode (see left panel of Fig. [2] (b)). 
In addition, at higher GR a significant Raman intensity is observed between the G and D 
modes, which originates from the presence of disordered carbon bondsP^Ell At lower GR 
(below the dashed line of Fig. |2| the peak intensity ratio 1(D) /1(G) greatly diminishes and 
the G and D bands become better resolved due to the Raman intensity between those two 
modes decreasing drastically. Both observations point to a larger crystal size and a higher 
crystal quality for lower growth rates. 

The two Raman spectra (GR=1.08 A/min), shown at the bottom of Fig. [2](b), reveal that 
by changing the excitation laser wavelength from A = 532 nm (green trace) to A = 633 nm 
(red trace) there is a clear redshift in the positions of the D and 2D Raman bands. In 
crystalline graphene layers, such a redshift arises from the wave-vector dispersion of the 
optical phononsPH^The size of the frequency shift that we observe is comparable to those 
reported for graphite and graphene P ^ 24 * 28 * The observed energy dispersion of the MBG films 
provides further evidence of crystallinityP^ 
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2. NEXAFS spectroscopy 



The two growth regions have the distinct NEXAFS signatures, as shown in Fig. |2](c). In 
both regions we find spectral fingerprints of sp 2 -hybridized carbonj-^ strong peaks at 285.4 
and 292.0 eV that correspond to excitation of a carbon Is core electron to the unoccupied 
7r* and a* orbitals, respectively. The sharpness of the NEXAFS features indicates a well- 
defined bonding environment and long-range periodic order in the electronic structure. The 
cr* fine structure, in particular, is specifically characteristic of graphite, and includes a sharp 
onset due to an excitonic core hole- valence state interaction and the broader a* peak at ~1 
eV higher photon energy due to more delocalized a* states.^ Thus, the NEXAFS spectra 
unequivocally prove the formation of sp 2 bonds between carbon atoms in the MBG films. 

NEXAFS is also sensitive to substrate-relative bond-orientations.^ Being governed by the 
transition dipole matrix element between a core electron and an unoccupied orbital above 
the Fermi level, the NEXAFS intensity depends upon the angle between the electric field 
vector of the incoming x-ray beam and the molecular orbitals in the system (see inset of Fig. 
[3]) . Hence, we directly probe the degree of bond anisotropy in the sp 2 films by changing the 
angle of the incident x-ray beam from near parallel (20°) to near perpendicular (70°) to the 
substrate; the E-field vector is perpendicular to the beam axis. 

For higher GR (upper half of Fig. [2] (c)) no angular dependence of the NEXAFS res- 
onances is observed, indicating a fairly isotropic arrangement of sp 2 bonds. In striking 
contrast, the NEXAFS intensity becomes strongly dependent on incident angle at lower GR 
(lower half of Fig. [2] (c) ) . The intensity of the n* (cr*) peak is at its maximum (minimum) 
at 20° and minimum (maximum) at 70° incidence, indicating highly oriented planar C=C 
bonds parallel to the substrate surface. Here, the sp 2 carbon layers grow in a two-dimensional 
plane. The ability to grow sp 2 carbon layers well aligned to the plane of the substrate, and 
the presence of two regions with distinctly different degrees of bond anisotropy is emphasized 
by the inset of Fig. [2] (c), which plots the area of the 7r* peak as a function of the incident 
angle for the two regions. 

The homogeneity of the material throughout the volume can be probed with NEXAFS 
by varying a bias voltage applied to the sample. By changing the voltage from -250 to 
-50 V, the depth within the carbon film from which detected electrons are emitted can 
be tuned from about 1 nm to about 7 nm (maximum film thickness 0o < 3.5 nm). The 
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higher voltage allows detection of electrons only from the near surface-region. Since we 
do not observe a bias-dependency of the spectral features, we conclude that the films are 
homogenous throughout the volume. This excludes the possibility of initial formation of a 
planar film in the isotropic region of the films followed by accumulation of defects as the 
film thickness is increased. 

B. Angle-dependent NEXAFS: results on various substrates 

Figure [3] compares angle-dependent NEXAFS spectra for few-layer (~2 nm) MBG 
graphene on SiC>2 and on mica, with a single high-quality graphene layer grown on Cu 
foil by CVD. 30 The similarity of the NEXAFS spectra from the MBG and CVD grown 
graphene is striking. In fact, aside from a slightly weaker angular dependence of the MBG 
films, the main difference between the MBG and CVD spectra is the intensity in-between 
the 7r* and the a* resonances, which is due to C-0 and C-H bonds (a resonance due to an 
interlayer state in few- layer graphene also appears in this region) EMS] 

This intensity between the n* and the a* resonances can be explained by the larger 
number of dangling bonds available at the grain boundary of the MBG nanocrystals, due to 
their smaller grain size compared to those in the CVD samples. These are readily saturated 
by oxygen and hydrogen bonds. These bonds tend to distort the planarity of graphene films 
and therefore may also explain the suppressed angular dependence of the NEXAFS data for 
the MBG films compared to CVD graphene. 

Importantly, no features associated with sp 3 carbon-carbon bonds are observed in the 
NEXAFS data. 

C. NEXAFS results from thick MBG films 

To further investigate the impact of disorder, we compare NEXAFS data from a thick 
MBG film (Go = 54.4 nm), both with that from glassy carbon (used as carbon source) and 
from highly-ordered pyrolytic graphite (HOPG). This is shown in Fig. |lj 

While the NEXAFS spectrum of the MBG film is very similar to that of HOPG, distinct 
differences are observed from the glassy-carbon spectrum, which has significant sp 3 content. 
In particular, the sp 2 — n* and o* peaks are strongly suppressed and the sp 2 — a* peak is 
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significantly broadened. An onset and a peak centered around 289 eV appears due to the 
sp 3 — cr* absorption edge of diamond and a C-H resonance J 32 l 33 l 

In contrast, the thick MBG film and the HOPG traces possess the spectral signatures 
of sp 2 bonds, as already discussed. HOPG has better long range periodic ordering, as is 
evidenced by the sharpness of the a* resonance. As in Fig. |3j the MBG films show some 
C-H and C-0 bonds at the grain boundaries of the nanocrystals, as well as non-uniform 
bonding between the differently oriented graphene nanocrystals in three dimensions, giving 
rise to the intensity between the sp 2 — a* and the it* resonances (indicated by the arrow in 
Fig. §. 

Thus, the thick MBG films are indeed structurally very similar to HOPG and consist of 
graphite nanocrystals. 

D. Scanning tunneling microscopy 

As expected for sp 2 bonded carbon, the MBG films show electrical conductivity at room 
temperature. Preliminary 4-probe transport measurements reveal a sheet resistivity of a 
few kO; sufficient conductivity for STM measurements. Figure [5] (a) shows a 3-dimensional 
ambient STM topography of a MBG film on a mica substrate (the size of the image is 
4x4nm 2 ). Several flat terraces are observed. A line profile, along the blue line in [5] (a), 
reveals 0.33 nm high steps, as shown in Fig. [5] (b). These step-heights are comparable 
to the interlayer distance in graphite,^ as one would expect in graphene multilayers. The 
surface roughness is dominated by the roughness of the underlying substrate. This has been 
confirmed by tapping-mode AFM measurements. 

IV. ANALYSIS OF RAMAN SPECTRA: NANOCRYSTAL SIZE 

Detailed analysis of Raman lineshapes enables estimates of the crystallite sizes. Typical 
Raman spectra of MBG nanocrystals grown on Si02 are shown in Fig. |6} A (linear) 
background subtraction between 1900 cm -1 and 2300 cm -1 , has been applied. 

Figure [6] (a)- (d) show results at high GRs. The blue shift of the G band to ~ 1600 cm -1 
(from ~ 1585 cm -1 in graphite) seen in Figs. [6] (a) and (c) is attributed to the unresolved 
superposition of the G and D' Raman modes. EHMH] rp^g j^ aman intensity between the D 
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and G band is tentatively interpreted as from disorder at the grain boundaries. In addition 
to the 2D (G') band at (~ 2672 cm -1 ) we also observe two second-order bands at higher 
Raman shifts (Figs. ||(b) and (d)): G+D at ~ 2928 cm" 1 and G+D' at ~ 3202 cm" 1 . These 
allow a calculation of the energy shifts of the G and D' modes, which thereby are found to 
be at 1584 cm -1 and 1618 cm -1 , respectively. All these values are in very good agreement 
with previous reports for such graphene-like systemsP^EH 

Figure [6] (e) shows results at lower GRs: Here the G mode redshifts to ~1585 cm -1 , 
indicating that the contribution of the D' band is reduced, most likely as a consequence of a 
larger nanocrystal size. In addition, the D-mode intensity is reduced relative to the G mode 
and the intensity between those two modes decreases. Four Lorentzians, corresponding to 
the three graphene optical phonon frequencies: D, G and D'; and a fourth one related to 
the 3TO Si phonon (at 1450 cm -1 ), reproduce the dataP^ 

The intensity ratio 1(D) /1(G), provides an estimate of the crystallite dimensions) 19 * 20 * 23 * 24 * 36 ! 
The graph of Figure [6] (f ) reveals an unambiguous trend: that the grain size increases up 
to 22 nm on reducing the GR. This result is consistent with the reduced Raman intensity 
between the D and G lines in Fig. [6] (e). 

V. CONCLUSION 

We have investigated non-epitaxial growth of graphene on insulating substrates by using 
a molecular beam of carbon atoms and have explored the influence of the key parameters 
required to obtain good quality, ultra-thin films by MBG. 

Using NEXAFS and Raman spectroscopy we have shown that growth rate is a crucial 
parameter for two-dimensional (layered) growth of graphene crystals, as it strongly influences 
the alignment of the sp 2 -bonds. NEXAFS spectra for high growth rates reveal isotropic 
orientation of the sp 2 -bonds. This growth can be regarded as quasi-3D. 

On reducing the growth rate we concomitantly increase crystallite size to ~ 22 nm and 
align the graphene multilayer-crystals parallel to the substrate. The reduction of grain 
boundaries manifests as reduced Raman scattering intensity between the D and G bands 
and anisotropy in the bond-orientations in angle- dependent NEXAFS measurements. 

Typical MBG parameters (growth rate, substrate temperature, surface mobility), and the 
MBG setup itself, offer a wide parameter space in which to explore the growth of a range 
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of layered materials with van-der Waals coupling between the layers. At the same time 
MBG allows for the potential growth of heterostructures based on these layered materials. 
We expect that even the use of smoother and more inert substrates, like hexagonal boron 
nitride, could greatly improve the crystal quality. 
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FIG. 1. (a) Schematic of the growth arrangement - the solid carbon source is placed closely facing 
the substrate. This geometry leads to a gradient in growth rate along the long axis of the sample, 
(b) Photograph of an ultra-thin graphene MBG film on a SiC>2 substrate. The changes in color are 
from the graded thickness, and the contrasting colour at the top is due to mechanical clamping of 
the substrate during growth, (c) Schematic of the MBE-inspired growth. 
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FIG. 2. (a) Schematic separating two distinctively different regions: High growth rates (HGR) 
result in isotropically arranged graphene layers. Low growth rates (LGR) are dominated by two- 
dimensional growth parallel to the substrate surface, (b) Color map of the Raman intensities along 
the GR gradient and concomitant varying growth rate. Typical features observed can be linked 
to graphitic material, such as the D, G, 2D(G') and G+D peaks. Inserted are two Raman traces 
at two excitation wavelengths at GR = 1.08 A/min. (c) NEXAFS spectra showing characteristic 
sp 2 bond features: ir* and a*. For HGR, no polarization dependence of the it* and a* resonances 
is observable. For LGR, a large polarization dependence is observed indicating layered growth of 
the graphene crystals. The inset depicts the dependence of the integrated ir* peak intensity on the 
beam orientation for both HGR and LGR. 
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FIG. 3. Angle dependent NEXAFS measurements in the low-growth-rate region. The similarity of 
MBG graphene grown on SiO"2 and mica with CVD graphene, including the energy position and 
angular dependence of the NEXAFS features, is clearly observable. Incident angles of 20° and 
70° correspond roughly to out-of-plane and in-plane polarizations, as shown schematically in the 
inset. 
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FIG. 4. Orientation independent NEXAFS of bulk material measured near the magic angle (50°). 
The spectra from HOPG and a 50 nm-thick layer grown by MBG show the same detailed structure. 
In contrast to that, the glassy carbon trace shows a strongly suppressed sp 2 -7r* resonance, a shift 
and broadening of the sp 2 -cr* and a well pronounced peak at ~289 eV, which is associated with 
the a* resonance of sp 3 -carbon. 
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FIG. 5. Ambient STM measurements reveal terraces with a step height of roughly 0.33 nm as 
shown in (b), which shows the line profile along the blue line in (a). 
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FIG. 6. Micro Raman spectra of MBG graphene nanocrystals on amorphous Si02 measured at 
various GRs: (a and c) the D and G modes, (b and d) the G'(2D), G+D and G+D' modes. The 
spectra in (c) and (d) are typical for the crossover region, (e) a micro-Raman spectrum acquired 
at low GR: The D'-mode intensity decreases and the region between the D and G bands changes 
to show an additional peak attributed to a 3TO-Si phonon. (f) The crystal grain size estimated 
from the ratio of the D and the G modes. Reduction of the GR leads to an increase of grain size. 
Scale bars represent one arbitrary unit. 



20 



